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ABSTRACT 

The next generation of telescopes aim to directly observe the first generation of galax- 
ies that initiated the reionization process in our Universe. The Lya emission line is 
robustly predicted to be the most prominent intrinsic spectral feature of these galaxies, 
making it an ideal target to search for and study high redshift galaxies. Unfortunately 
the large Gunn-Peterson optical depth of the surrounding neutral intergalactic medium 
(IGM) is thought to render this line extremely difficult to detect prior to reionization. 
In this paper we demonstrate that the radiative transfer effects in the interstellar 
medium (ISM), which cause Lya flux to emerge from galaxies at frequencies where 
the Gunn-Peterson optical depth is reduced, can substantially enhance the prospects 
for detection of the Lya line at high redshift. In particular, scattering off outflows of 
interstellar H I gas can modify the Lya spectral line shape such that >5% of the 
emitted Lya radiation is transmitted directly to the observer, even through a fully 
neutral IGM. It may therefore be possible to directly observe 'strong' Lya emission 
lines (EW >50 A rest frame) from the highest redshift galaxies that reside in the 
smallest H II 'bubbles' early in the reionization era with JWST. In addition, we show 
that outflows can boost the fraction of Lya radiation that is transmitted through the 
IGM during the latter stages of reionization, and even post-reionization. Coupled with 
the fact that the first generation of galaxies are thought to have very large intrinsic 
equivalent Lya equivalent widths, our results suggest that the search for galaxies in 
their redshifted Lya emission line can be competitive with the drop-out technique out 
to the highest redshifts that can be probed in the JWST era. 

Key words: galaxies: high redshift - (galaxies): intergalactic medium - ISM: kine- 
matics and dynamics ~ line: profiles - scattering ~ radiative transfer 
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1 INTRODUCTION likely top-heavy (|Larsonlll99llBromm et al.ll2002l '). Both the 

top-heavy IMF and low (or zero) gas metallicity enhanced 
One 01 the mam science drivers oi the next feneration , „. . . , , ,i , n , i • 

^ , . , , ^ . ? , . m the number oi ionizing photons that the nrst galaxies emit- 

oi telescopes is to detect the farst generation ot galaxieaj i.,i,ri i i - . ni. r 

ted compare d to that ot local galaxies , at a faxed star tor- 



that formed in our Universe. These galaxies contained hot- ^. ^ ir; c oi 7il UrTAnl 'it; — I — nUorvil 

, r-^ . ^, mation rate ( Tumhnson fc Shulll 12000 : IBromm et al.l 120011 : 

ter and more compact stars (e.g. I lumfanson & bhull 200(J; rT"; IITTT^ L^^oIn rni ■ \ ^ . . . . , 

nrr^^rrf , . . ■ / — ^ ^t^x — lachaerenliOO^ . iOOdl l. 1 his enhancement m lonizmg lummos- 

Bromm et al. 20011 ). whose initial mass function (IMF) was 



ity results in larger H II regions in the interstellar medium 
(ISM) of these galaxies. As a result, one of the key predicted 

* E-mail:mdijkstr@cfa.harvard.edu properties of the first galaxies are prominent nebular emis- 

1 The definition of 'the first generation of galaxies' is somewhat sion lines, dominated in flux by hydrogen Lya (A = 1216 A, 

arbitrary. We take it to mean any star forming galaxy during see e.g. Johnson et al. 2009). The first generation of galax- 

the earliest stages of the reionization that might be detected (in igg are therefore likely to have been strong Lya emitters, 

either their continuum or in on of their lines) by the next genera- with equivalent widths possibly as high as EW~ 1500 A 

tion of telescopes. Although the analysis presented in this paper (|Schaeredl200l [jooj : Ijohnson et akllioogl . also see Partridge 

also applies to the first stars that likely formed one-by-one in P p Ki i q«7 i\/r ■ i Q7R\ 

mimnalos, the overall fluxes from these sources is much too taint ' 

to be detected in the near future (regardless of radiative trans- In this paper we investigate the prospects for detect- 

fer effects that may boost the detectabihty of Lya emission from ing this Lya emission. The first galaxies were surrounded 

such sources). by a mostly neutral intergalactic medium (IGM), which is 
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extrem ely optically thick to Lya radiation. iLoeb fc Rybickil 
l|l999l) showed that scattering of Lya photons through a neu- 
tral IGM caus es galaxies to be surrounded by diffuse Lyg 
halos (also see iKobayashi fc Kamayall2004l : iKobavashi et al.l 
l2006l ). While the total flux in these halos can be substantial, 
their large angular size results in surface brightness levels 
that are beyond reach of even future telescopes such as the 
James Webb Space Telescope (JWST, see §[2TTJ- 

If Lya radiation from the first galaxies is to be detected, 
it will therefore be via Lya photons that were transmitted 
directly to the observer. In a neutral IGM, there are two 
mechanisms by which this can be achieved: 
(i) The first mechanism is due to ionizing radiation that 
escapes from the galaxy, which creates a surrounding H II 
'bubble' that can strongly boost the detecta bility of Lya 
emissio 

rE (e.g. lHaimanll2003 : ICen et ahllioosl ). This is be- 
cause Hubble expansion redshifts Lya photons while they 
propagate freely through the ionized gas. As as result, a 
fraction of photons enter the neutral intergalactic medium 
on the red side of the line center, where the IGM optical 
depth optical depth can be smaller by orders of magnitude. 
For example, th e Gunn-Peterson optica l depth at redshift z 
IS given by (e.g. iBarkana fc LoeblbOOlh 
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where xm denotes the neutral volume fraction of hydrogen 
in the I GM. The Gunn-Peterson optical depth (Eg [1]) re- 
duces to (|Miralda-Escuddll998l : [biikstra fc Wvithell2006l ) 
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^600 kms~V V 10 
for photons that enter the neutral IGM with a redshift of Av 
from line center. Eq [2] illustrates that photons that enter the 
neutral IGM significantly redward of the Lya resonance are 
scattered with only a weak optical depth. However given the 
high density of the surrounding IGM it is generally assumed 
that the first galaxies would not have had large enough H II 
'bubbles' to prevent complete damping of the line. 
(ii) The second mechanism for enhancing the transmission 
of Lya photons from high redshift galaxies is due to radia- 
tive transfer effects within the ISM of galaxies (in particular 
outflows of interstellar H I gas), which can shift Lya pho- 
tons to the red side of the line before it reaches the IGM. 
Observed interstellar metal absorption lines (Si II, O I, C II, 
Fe II and Al II) in Lyman Break Galaxies (LBGs) are typ- 
ically strongly redshifted relative to the galaxies' systemic 
velocity (with a median off-set of ~ 160 km s~^), while the 
Lya emission line is strongly re dshifted (with a m edian ve- 
locity offset of ~ 450 km s~^ Isteidel et al.ll201ol . also see 
Shapley et al. 2003). Thi s suggest that large scale outflows 
are u biquitous in LBGs (|Shapley et al.l l2003l : ISteidel et al.l 
I2OIOI ). Furthermore, scattering of Lya photons by H I in 



^ Indeed, the fact that the reionization process likely affects 
the observed number and distribution of high-redshift Lya 
emitting galaxies is exactly the reason why Lya emitting 
galax ies are thought to pr ob e this e poch (Haiman S 
19991: iMalhotra fc Rhoadd |2004 iFurlanctto ct al 
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outflows can successfully explain ob served Lyg line shapes 
in Lyg emitting galaxies a t 2 = 3 — 6 jVerhamme et al" ] |2006l . 
|2004 IVanzella et al.|[2009h . 

In this paper we explore the outflow mechanism in more 
detail. We will show that these 'local' (i.e. inherent to the 
galaxy itself) processes in the ISM can cause as much as 
>5% of the emitted Lya radiation to be directly transmitted 
to the observer even through a fully neutral IGM. This result 
is important for the study of high redshift galaxies, because 
it suggest that detecting the Lya emission line of the first 
generation of galaxies may well be within reach of the next 
generation of telescopes including JWST. 

The outline of this paper is as follows: we describe our 
models and present our results in § (2] We discuss our models 
and the implications of this work in § |31 before we conclude 
in § m The cosmological parameter values used throughout 
our discussion are (n ^, ^a, ^i,, h) = (0.27, 0.73, 0.046, 0.70) 
l|Komatsu et al.ll2009l ). 



2 MODELLING OF LYg FROM HIGH 
REDSHIFT GALAXIES 

We divide our study into modelling of three phenomena 
which are important in the apparent brightness of Lya emit- 
ters at high redshift. Firstly in § 12.11 we compute the in- 
tegrated (over the sky) flux density, and integrated (over 
frequency) surface brightness profiles of the Loeb-Rybicki 
halos that surround the first first generation of galaxies. 
This calculation illustrates the difhculty of directly detecting 
these Lyg halos. Then in ij |2.2l we show how radiative trans- 
fer through a static ISM causes the predicted Loeb-Rybicki 
halos to be accompanied by Lya point sources, which are 
more easily detectable. Finally in § 12.31 we show that this 
result is strengthened further when the ISM of galaxies con- 
tain outflows of H I. In all our calculations we assume that 
the intrinsic, i.e. the emitted, luminosity of the galaxy is 
La = 10*^ erg s~^, which corresponds to a st ar forma- 
tion rate of SFR <,2Mq yr"^ for Z ^10"^Zq (|Schaerej 
|2003| . for a Salpete'r IMF from 1 - WOA'Iq). For compari- 
son, the brightest observed Lyg sources at z — 6.5 — 6.6 
have 'observed' luminosities (i.e. the observed flux times 
4-Kdj;^{z), where dL(z) denotes the lu minosity distance to red- 
shift z) of a few ti mes 10'*'^ erg s~^ l|Kashikawa et al" I I2OO6I : 
lOuchi et al. II2OI0I ) . Depending on what fraction of the emit- 
ted Lya flux we detect, their intrinsic luminosities can be 
significantly higher. Our results scale linearly with La- 

Throughout this paper, we denote photon frequency 
with the dimensionless parameter x — — Va)IAva, 
in which Va = 2.47 x 10^"" Hz denotes the Lyg rest- 
frame frequency; Ava = i^a^th/c, fth = \/2kBT /nip = 
12.9(r/10'' K)^/^ km s"\ where T denotes the gas temper- 
ature. Here, c, ks, and nip denote the speed of light, Boltz- 
mann constant, and the proton mass, respectively. In all our 
calculations we assume that the gas temperature in the ISM 
is TisM = 10* K, but we have explicitly verified that our fi- 
nal results do not depend on this assumption. Furthermore, 
the Gunn-Peterson optical depth for photons that enter the 
IGM while in the red wing of the Lya line is independent of 
temperature (see Eq[5)), and the assumed IGM temperature 
of TiGM = 300 K (appropriate for neutral intergalactic gas, 
see e.g. Pritchard & Loeb 2008) is irrelevant. 
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Figure 1. The observed properties of Lya halos (a.k.a Loeb-Rybicki halos) that surround galaxies embedded within a fully neutral 
comoving IGM. The 'noisy' features in these plots are due to the finite (here Afp^ot = 10^) number of photons used in our Monte-Carlo 
calculations. Left panel: the histogram shows the integrated (over the entire area on the sky) emerging spectrum, under the assumption 
that all Lya photons were emitted at line center, and assuming an emitted Lyo luminosity of La = lO'*'' org s~^. The lower horizontal 
axis shows the dimensionless frequency parameter x, while the upper horizontal axis shows the frequency parameter that was employed 
by Loeb & Rybicki (1999). The peak flux density occurs at x ~ —300 which in our model corresponds to a redshift of 660 km s~^. 
The FWHM of the spectrum is ~ 1500 km s~^. The vertical axes are expressed in arbitrary units on the left, and physical units on the 
right. The peak integrated flux density reaches 30 njy. For comparison NIRSpec aboard JWST is expected to reach a sensitivity of ~ 
hundreds of nJy in 10* s (S/N=10) for R=2700, or line fluxes of ;^10-l^ erg s ^ cm ^. Right panel: the integrated (over frequency) 
surface brightness 5 as a function of impact parameter (9) from the galaxy. We find that S < 10~^^ erg s~^ cm"'^ at all 9. This plot 
illustrates that it will be extremely difficult to directly detect these Lya halos (see text). 



2.1 Loeb-Rybicki Halos 

Figure [T] shows the observed properties of Lya halos that 
surround galaxies embedded within a fully neutral comov- 
ing IGM. The histogram in the left panel shows the inte- 
grated (over the entire area on the sky) emerging spectrum 
computed wi th the Monte-Carlo L ya radiative transfer code 
McHammer (|Diikstra et al.ll2006D . Following previous work, 
we assumed in this calculation that the Lya photons were 
emitted at line center, i.e. a; = 0. The lower horizontal axis 
shows the dimensionless frequency parameter x, while the 
upper horizontal axis shows the frequency parameter that 
was employed by Loeb & Rybicki (1999). All Lya photons 
emerge with a systemic redshift (i.e. a; <C 0). The peak flux 
density occurs at x ~ —300 which corresponds to a red- 
shift of ~ 660 km s"\ The Full Width at Half Maximum 
(FWHM) of the spectrum is ~ 1500 km s"^ The vertical 
axes contain arbitrary units on the left, and physical units 
on the right. The peak integrated flux density is ~30 nJy. 

The right panel shows the integrated (over frequency) 
surface brightness S (in erg s"'^ cm~^ arcsec"^) as a func- 
tion of impact parameter {8, angular separation) from the 
galaxy in arcsec. The surface brightness contains a core out 
to S ~ 20 arcsec after which it drops steadily. We find 
that S < 10~'^^ e rg s~^ cm~^ at all 6. For comparison, 
I Ranch et al.l (|2008l ) reached a 1— a surface brightness limit of 
•S'lim = 8 X 10"'^" erg s~^ cm~^ in a 92 hr long exposure with 
the ESO VLT-FORS2 instrument, which represent the deep- 
est observations to date. This, when combined with JWST's 
proposed sensitivity limil|f| to point sources of ^10~^^ erg 
s~^ cm~^, implies that it is and will remain extremely diffi- 
cult to detect these Lya halos, even with the next generation 
of telescopes. 

^ See http://www.stsci.edu/jwst/scionce/sensitivity/ This esti- 
mate assumes R = 1000 or ij = 2700, and assumes an integration 
time of 10^ s. 
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2.2 Static ISM 

In the previous calculation we assumed that all Lya escaped 
from the galaxy at the line center. In reality complex radia- 
tive transfer effects occur inside the ISM of a galaxy. H I ob- 
servations show that local galaxies contain H I column den- 
sities of A^Hi ~ 10^*^ - 10^^ cm . Such large columns of H I 
gas can affect the spectrum of Lya photons as they emerge 
from the galaxy. Furthermore, we expect column densities 
of gas in high redshift galaxies to increase as oc {1 + z)^ 
(see the end of this section for a more quantitative discus- 
sion). It is therefore plausible that Lya photons need to 
traverse a substantial, i.e. A''hi > 10^" cm~^, column of H I 
gas before escaping from the galaxy. To estimate the impact 
of large H I column densities on the detectability of Lya 
emission, we turn to analytic solutions for the spectrum of 
Lya photons emerging from static, homogeneou s, extremely 
optica lly t hick med i a. FoU o wing the analy s es of iHarringtonl 
( 1973[ ) and lNeufeldl l|l990^ ■ lDiikstra et al.l (|2006l ) computed 
the Lya spectrum emerging from a sphere to be 




assuming a source of photons in the center of the sphere. 
Here, tq = 5.9 x 10'^ {Nhi/W'^° cm-^){TisM/10^ K^^^ de- 
notes the fine center optical depth from the center to the 
edge of the sphere, and a = 4.7 x 10"*(risM/10'* K)"^/^ 
denotes the Voigt parameter. The function J(x) is nor- 
malized to (27r)~^, and its maximum occurs at Xp ± 
0.9(a„ro)i/3 ~ ±13{Nm/W^° cm-^y^^TisM/W^ K)'^^\ 
This dimensionless frequency shift translates to Av = 
±160(ArHi/10^° cm-2)^''3(risM/10* K)^'"^ km s-\ The 
weak temperature dependence (oc T^^^) explains the weak 
dependence of our results on the assumed ISM temperature. 
The Lya spectrum emerging from a semi-inflnite slab with 
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Figure 2. The solid histogram in Figure shows the spectrum 
of Lya photons that were not scattered at all. These unscattered 
photons would be observed as a point source. The grey line shows 
the intrinsic flux density of the source, Jint(a;), which corresponds 
to the observed Lya spectrum if we were able to detect all Lya 
radiation (EqlSjl. Overplotted as the black solid line shows the 
analytic calculation of the observed spectrum JpointC^^) given by 
Jint(x) exp(— TjQM [2^]) (see text). The inset shows the fraction of 
directly transmitted Lyo flux, /trans, as a function of H I column 
density A^hi. The transmitted fraction increases from ~ 1% at 
ATjjj = lo2i cm-2 to ^5% at A^hi ^7 x lO^^ cm'^. A trans- 
mission /trans of Only a few per cent is important: if the intrinsic 
Lya EW of the first generation of galaxies was as high as 1500 A, 
then the observed restframe EW is EW~ 50(/trans/0.03) A(see 
text) 

the same line center optical depth is broader by a factor of 
1.1, but is otherwise very similar. 

We compute how the predicted properties of the ob- 
served Lya radiation change when we inject photons into 
a neutral IGM following the distribution given by Eq [31 
rather than in the line center (as in § I2.1|l . The most im- 
portant difference is that a non-zero fraction of the photons 
are transmitted directly to the observer without scattering). 
The solid histogram in Figure [2] shows the spectrum of Lya 
photons that were not scattered at all. These unscattered 
photons would be observed as a 'point source' (the angular 
scale is set by the physical scale of the scattering medium), 
and its peak flux density is ~ 4 njy. This peak flux den- 
sity is smaller than the peak integrated flux density of the 
Loeb-Rybicki halos, because in this model the total fraction 
the photons that was transmitted to the observer directly, 
/trans, was Only ~ 1%. For comparison, we have overplotted 
the intrinsic flux density of the source as the grey line, as 

* The scattered radiation is spread out in a halo that resembles 
the Loeb-Rybicki halos. However, at small angular separations 
(5 gl arcsec) we find that the surface brightness profile is en- 
hanced by a factor of ~a few-ten (see Appendix |A] for an exam- 
ple). This is still well below detection thresholds of existing and 
future facilities. Furthermore, this enhancement in the central sur- 
face brightness profile vanishes when we allow for the existence 
of an H II region around the galaxy (see Fig lAll l. 



given by EqO This corresponds to the observed Lya spec- 
trum if we were able to detect all Lya radiation. Overplotted 
{black solid line) is the analytic calculation of the observed 
spectrum Jpoint {x) given by 

Jpoint(a;) = J(a;) exp(-riGM[a;]), (4) 

where J{x) is given by Eq[3] Furthermore, riGM(a;) denotes 
the opacity of the neutral IGM to a Lya photon that es- 
capes from the galaxy at frequency x, whic h is riGM(a^) = 
TGP.o-^ J^^4>{x')dx' (e.g. Appendix A of iDiikstra et all 

Here, tgp.o denotes the Gunn- Peterson optical 
depth given by Eq [T] and <t>(x ) denotes the Voigt func- 
tion (|Rvbicki fc LightmanI [l979l ) . Not suprisingly, the an- 
alytic calculation agrees well with the Monte-Carlo calcula- 
tion (which in the end should generate random realizations 
of this function). We point out that the predicted Lya spec- 
tral line shape is not as asymmetric as is often observed 
in lower redhsift galaxies. This is because the Lya absorp- 
tion cross-section varies weakly with frequency (as Ad~^, 
see Eq[2]) in its damping wing, and scattering in the IGM 
introduces no 'sharp' cut-off in the observed Lya line shape. 

The mset of Figure [2] shows how the fraction of directly 
transmitted Lya flux depends on H I column density A^hi. 
As mentioned previously, /trans ~ 10~^ at A^hi = lO^'^ cm~^. 
The transmitted fraction increases to ^5% at A^hi <;7 x 
lO'^^ cm~'^. For comparison, the central column density for a 
standard Mo-Mao & White exponential disk-when seen face 
on-is Nnifi ~ 3 x 10^^([1 + z]/ll)^/^(ucirc/13 km/s), under 
standard assumptions that the spin parameter A = 0.05, 
and that jd = md = 0.05. This suggests that H I column 
densities in excess of lO'^^ cm~^ are expected at the redshifts 
of interest, and that therefore that /trans can exceed a few 
per cent. 

We stress that although small, a transmission /trans 
of only a few per cent is important. The total observed 
flux in the Lya point source at 2 = 10 is faha = 3.7 x 
10"^" (/trans/0.05) ~ 10"^^ erg s"^ cm-2 which is within 
reach of JWST. For example, an R = 1000 grating would 
yield a 5 — (T(/trans/0.05) detection in a 10^ s exposure 
(band I). In addition to rendering the line emission from 
very high redshift galaxies detectable by JWST, this effect 
could lead to strong Lya, since given an intrinsic Lya EW 
for the first generation of galaxies of EWint = 1500 A, the 
observed restframe EW is EW~ 50(/trans/0.03) A. For com- 
parison, only ^4% of LBGs at z = 3 have emission lines 
with EW ^50 A (jShaplev et al.ll2003l ). Indeed, for an ob- 
served rest frame EW~ 75(/trans/0.05) A, and a line flux 
of /obs = 3.7 X 10-1'' (/trans/0.05) ~ lO"!'' erg s-i cm-2 

^ This implicitly assumes that the gas in IGM is at mean density 
of the Universe, and that it is co-moving with the Hubble flow. In 
reality we expect strong departures from the Hubble flo w in close 
proxi mity to halos, and that the gas is overdense there l lBarkanal 
|2004) . If we incorporate this into our model, then this increases 
the opacity of the IGM. However, ionizing radiation that escapes 
from the galaxy ionizes most of this denser gas. Indeed, we have 
verified that for escape fractions of ionizing radi ation as low as 
.fesc = 1%, more realistic models of the IGM (as in lDiikstra et all 
l2007bl ) actually transmit more Lya. The calculations that are 
present in this paper are therefore conservative. 
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Figure 3. The black solid line shows the spectrum of Lya pho- 
tons that were transmitted directly to an observer through a fully 
neutral IGM. The grey solid line shows the intrinsic spectrum 
of this galaxy (see Ahn et al. 2003 and Verhamme et al. 2006 
for a detailed discussion on these features in the spectrum). The 
ISM of this galaxy was modeled as a thin, outflowing (with speed 
^'sh)i spherically symmetric shell of H 1 gas (with column den- 
sity A^Hi)- This type of model successfully reproduces observed 
Lya line shapes in known Lya emitting galaxies. We assumed 
that (Afni.-Ush) = (10^'^ cm-^, 200 km s"!). Because a signifi- 
cant fraction of the radiation comes out of the galaxy with a 
large systematic redshift, /trans ~ 4% of all emitted photons is 
transmitted directly to the observer. 
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Figure 4. This figure shows the fraction of Lya photons, /trans, 
that is transmitted directly through a fully neutral IGM at 2 = 10, 
as a function of v^i-^ for A'hi = 10^" cm~^ (red squares) and 
A^Hi = 10^^ cm~^ (blue squares). Outflows of ^50 km can 
boost /trans tremendously to values as large as /trans ~ 10 — 20% 
(for Vah = 200 km s'^). That is, 10 - 20% of all Lyo photons 
may be transmitted through a fully neutral IGM. 

(typical for lower redshift Lya emitters), the corresponding 
continuum flux density is fi, ~ 2.7 njy, which would only 
yield a ~ 2 - cr detection in lO'^ s with NIRCAM on JWST. 
The Lya line flux may therefore be (slightly) more easily 
detectable than the continuum emission, even when we only 
observe a few per cent of the emitted Lya flux. 
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Figure 5. Same as Figure [J] but at 2 = 


15. 



2.3 Outflowing ISM 

Observed galaxies typically show evidence that large scale 
outflows are ubiquitous in LBGs: interstellar metal absorp- 
tion lines (Si II, O I, C II, Fe II and Al II) in Lyman 
Break Galaxies (LBGs) are typically strongly redshifted rel- 
ative to the galaxies' systemic velocity (with a median off-set 
of ~ 160 km s~^), while the Lya emission line is strongly 
redshifted (with a median velocity offset of ~ 450 km s~^ 
ISteidel et ahlboiol . also see Shapley et al. 2003). Scattering 
of Lya photons by H I in these outflows has successfully 
explained observed Lya line shapes -and their redshifts-in 
Lya emitting galaxies at z = 3 — 6 (| Verhamme et alll2006l . 
l2008l : IVanzella et al.|[2009l . see 5 13.51 for a discussion on this 
model) . 

We therefore repeat the exercise of § 12.21 for a suite 
of outflow models. Following ^Verhamme et al. (2006, 200^), 
we model the outflow as a spherically symmetric thin shell 
of gas that con tains an H I column density A^hi (also see 
lAhnet al.|[20o3 ). and outflow velocity Vsh- We assume that 
the shell has a radius of 1 kpc and a thickness of 0.1 kpc, 
but stress that the precise physical scale of the outflow is not 
important for our results. Our assumed gas temperature in 
the outflowing H I shell of Tism = 10* K corresponds to a b- 
parameter of 6 ~ 13 km s~^ in the terminology of Verhamme 
et al. (2008). Our results do not depend on this choice for 6, 
as the amount of flux at large Av depends very weakly on 
this parameter (see the right panel of Fig 15 of Verhamme 
et al. 2006). We further assume the H I shells to be dust- 
free (see § 13. 3p . IVerhamme et all (|2008l ) typically found that 
log Ahi ~ 19 — 22, and ~ — 500 km s~^, and this is 
the range of parameter space we explore. We compute Lya 
spectra emerg ing from the outflow s with the Monte-Carlo 
transfer code (|Diikstra et al.|[2006l ). In our calculations, the 
Lya photons are emitted at line center, i.e a; = 0. We veri- 
fied that our results are insensitive to this assumptiorjfl. We 
compute the impact of the IGM on the directly observed 



° We repeated some calculations in which Lya photons were 
emitted following a Gaussian distribution with a standard de- 
viation of (7 = 100 km s~^. While the precise emerging line 
shapes changed, the overall transmitted fraction remained prac- 
tically identical. 
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fraction of Lya by simply suppressing the intrinsic spectrum 
by exp(-TiGM[3;]) (see §[2]2}. 

The grey solid line in Figure [3] shows an example of the 
intrinsic spectrum for a model in which we assumed that 
{Nm,Vsh) = (10^° cm"^,200 km s"^). The intrinsic spec- 
trum is highly asymmetric, with more flux coming out on 
the red side of the Lya line center. The spectrum peaks at 
about ~ 2-i;sh, as expected for radiation that scatters back 
to the observer on the far side of the galaxy (see Verhamme 
et al. 2006, and Ahn et al. 2003 for a detailed discussion 
on these features in the spectrum). However, a significant 
fraction of the radiation comes out at larger redshifts. This 
radiation can be transmitted directly to an observer through 
a neutral IGM, and would be observed as a point-source. The 
black solid line shows the spectrum of the point source (see 
Appendixl^for a plot of the surface brightness profile of the 
scattered radiation). We find that the spectrum of the point 
source contains /trans ~ 4% of all emitted photons. Note 
that the FWHM of the observed spectrum of the galaxy, as 
well as the offset arise purely from radiative transfer effects. 

The directly transmitted fraction /trans depends on 
both A'hi and Wsh, In Figure|4]we show /trans as a function of 
Wsh for TVhi ~ 10^" cm~^ {red squares) and A^hi = 10^^ cm~^ 
[blue squares). Having outfiows as low as ^50 km s~^ makes 
a significant difference in the predicted /trans. For example, 
we find /trans ^5% for Ush ^200 (50) km s"^ for logiVni = 20 
(logAHi = 21). We stress that the outflow parameters were 
chosen to lie the range required to explain line shapes of ob- 
served LAEs at z < 6. If the dark matter halos that host the 
highest redshift galaxies had lower masses (as expected in a 
hierarchical structure formation scenario), then outflows at 
higher redshift may have occured at lower speeds, because 
of the observed (at lower redshifts) scaling of outflow veloc- 
ity with circu lar velocity of the host dark matter halo (e.g. 
iMartinllioOsI ). On the other hand, higher redshift galaxies 
are expected to be more compact, and therefore we naively 
expect larger column densities of H I gas at higher redshifts 
(see § 12. 2p . in which case low outflow velocities can boost 
/trans trcmcudously. The outflow properties in the highest 
redshift galax;ies, and their impact on Lya propagation is 
clearly a topic that needs further investigation. 



3 DISCUSSION 

In this section we discuss a range of issues arising from our 
model. 

3.1 Detecting Lya Emission During and After the 
Epoch of Reionization 

The results of the previous section focused on the flrst gen- 
eration of galaxies that were surrounded by a neutral IGM. 
However our work applies more broadly. For example, radia- 
tive transfer effects in the ISM can also boost the detectabil- 
ity of the Lya line emitted by galaxies during the later 
stages, and even post-reionization. It is generally thought 
that ionized intergalactic gas is transparent to Lya radia- 
tion. More speciflcally, residual H I that exists inside the 
ionized IGM is assumed to only affect radiation frequen- 
cies that lie blueward of the Lya resonance. However, in 
the standard cosmological model (ionized) intergalactic gas 



in close proximity to galaxies is expected to be overdense 
(1-1-5 ~ 2 — 20, Barkana 2004). Furthermore, gravity causes 
this gas to flow towards galaxies. This inflowing, denser ion- 
ized gas is expected to scatter a significant fraction of the 
Lya photons out of the line of sight. Even the ionized IGM 
therefore transmits on average as lit tle as ~ 10 — 30% of the 
Lya to an observer (iPiikstra et al.ll2007bl : llliev et al]|2008l : 
IZheng et al.ll2010al : lD'aval et aLlbOld ). It is therefore unclear 
whether H II regions during the reionization process actu- 
ally transmit enough Lya flux to an observer, to render it 
detectable. 

However in the presence of H I outflows, the transmitted 
fraction of Lya radiation through ionized gas can increase 
dramatically. This is because the outflowing H I gas imparts 
a redshift to the Lya photons that is well in excess of ~ a few 
hundred km s~^ (see i; l2.3l and Fig|3]). On the other hand, the 
influence of the ionized IGM only extends out to At; ^Winfaii 
(|Barkanall2004l l . Here, the inflow velocity, Winfaii, is of order 
the circular velocity of the dark matter halo hosting the 

galaxy, Winfall ^ Vdrc = 80(Mhalo/10'° Mo)'/3(^ _^ 2]/ll)i/2 

km s~^ (e.g. iBarkana fc Loebl l200ll . where Mhaio denotes 
the dark matter halo mass) . The impact of the ionized IGM 
in close proximity to the galaxy is therefore signiflcantly 
weaker when H I outflows are important in shaping the Lya 
spectral line shape that emerges from galaxies. 

This implies that winds can have important conse- 
quences when interpreting the observed sudden drop in 
the Lya luminosity b e tween z = 6.5 and z = 5.7 
(IShimasa ku et al. I l2006l : iKashikawa et al.1 l2006l : lOta et al.l 
|2008, also see Ouchi et al. in prep). The most important 
aspect of this observation is that the rest-frame UV luminos- 
ity function of these same galaxies does not evolve between 
these redshift s (within the uncertai nties, Kashikawa et al. 
2006). Indeed. [Piikstra et al.l (|2007al ) have shown that these 
two observations combined translate to a reduction in the 
number of detected Lya photons from z — 6.5 by a factor of 

1.1 — 1.8 (95% CL) relative to z = 5.7. The simplest interpre- 
tation of this observation is that the IGM at z = 6.5 is more 
opaque to Lya photons than at z = 5.7, be cause the IGM 
naturally only affects the observed Lya flux. iDiikstra et al.l 
(|2007ah argued that this can be explained quite naturally 
by an evolution in the opacity of the ionized IGM. However, 
if winds reduce the impact of the ionized component, then 
this conclusion becomes uncertain, and the observed reduc- 
tion in the Lya luminosity function may be at least partly a 
reionization-induced signature. In order to understand the 
evolution of the transmission in terms of the reionization 
history it is therefore crucial to understand the properties 
of winds in high-redshift galaxies, and how they affect the 
transport of Lya radiation. 

3.2 Detecting Other Lines 

In this paper we have focused our attention on the Lya 
emission line, because it is predicted to be the strongest 
spectral line of high redshift galaxies. It may be possible to 
detect other spectral lines as well. 

The most prominent alternative emission line is Ha 
(A = 6536 A), which is intrinsically weaker than Lya in flux 
by a factor of ~ 8, under the assumption of case-B recombi- 
nation. However only a fraction /trans = 0.05 of the emitted 
Lya radiation is transmitted, so that the observed Ha flux 
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may be larger by a factor of ~ 2.5(/trans/0.05)~^ . On the 
other hand the Ha Une lies deeper in the IR, A = 7([l+z]/ll) 
/xm, where JWST's Mid Infrared Instrument would yield a 
lower S/N detection of the Ha line for medium resolution 
spectroscopy. This suggests that it may be more difficult to 
detect Ha than naively expected. 

In addition the Helium Balmer a (HeH A =1640 A) 
line could contain a flux that is com parable to that of 
Ha l|Oh et aLlbOoH : I Johnson et al.ll2009l ). However, the pre- 
dicted strength of this line depends sensitively on the IMF, 
and it can be weaker than Ha by an order of magnitude for 
reasonable model assumptions l| Johnson et al.|[2009i V On the 
other hand, the composite spectrum of observ ed z = 3 LBGs 
contains a Hell A =1640 A emission line jShapl ev et al.l 
bOOj ). The Full Width at Half Maximum of this line is 
FWHM~ 1500 km s"\ which is broader than the FWHM 
of most other nebular emission lines. The origin of this line 
is not resolved. The li ne has been associated wi th popula- 
tion HI star formation (| Jimenez &: HaimanllioO^ ') occurring 
in pockets of pristine gas that persisted down to low red- 
shift, but may also originate in winds associated with mas- 
sive Wolf-Rayet stars. In this latte r case the line is expecte d 
to vanish for low gas metallicities l|Brinchmann et al]|2008h . 
Thus it is not clear whether the detection of Hell A =1640 
A in z=3 LBGs implies that this line should be present in 
the higher redshift galaxies. 

Finally, we expect numerous recombination lines asso- 
ciated with metals heavier than Helium (such as [O II], 
[O III], [S II], e.g. Huchra 1977). Again, robust pre- 
dictions of the strengths of these lines (which fall in the 
restframe optical) do not exist as it requires an accurate 
knowledge of the metallicity of H II regions in high-redshift 
galaxies. 



found that on average /d.csc ~ 10~^'^^'-^~^^ [Bouwens et al.l 
(2010) recently argued that candidate galaxies at z = 7 are 
considerably bluer than at lower redshifts. Their most lu- 
minous candidates are consistent with E(B-V)~ 0.05, which 
translates to /d,csc ~ 0.4 (which is significantly higher than 
Lya escape fraction inferred for galaxies at z <.2 see e.g. 
Hayes et al. 2007,2010). If we take this number literally, 
then our quoted Lya fiuxes are high by a factor of ~ 2.5 
(which implies that dust is not the dominant uncertainty). 
Moreover, in hierarchical models of galaxy formation, we ex- 
pect galaxies at higher redshifts, and lower luminosities to 
be less evolved, and hence the impact of dust is likely to be 

weaker. 

The relation derived bv lVerhamme et al.l l|2008l 'l was ob- 
tained by averaging over a dozen galaxies, and there ex- 
ists considerable scatter around their derived relation. For 
example, observations of low-redshift star forming galax- 
ies with the International Ultraviolet Explorer (lUE) and 
the Hubble Space Telescope (HST) revealed very weak Lya 
emiss ion (or even strong absorption) in metal po or ob- 
jects (|Hartmann et al.lll984 11983 : iKunth et al.lll994l '). while 
strong Lya emission co uld be detected from some metal and 
dust rich galaxies (e.g. iLeaueux et al]|l995f ). This scatter is 
likely related to the precise geometry of the dust distribu- 
tion JSc^^^^^Lj2009,), viewing angle towards the galax- 
ies (|Laursen et al.ll2009l ) and outflows. When outflows are 
present in galaxies, the strength of Lya emission appears to 
be independent of the dust content of the galaxy. In the ab- 
sence of strong outflows howe ver, Lya line strength appears 
to de crease with dust content (|Kunth et al ]|l998| -[ Atek et al.l 
l2008h . Our work implies that in high redshift galaxies, H I 
outflows further affect the subsequent impact of the IGM on 
the detectability of Lya emitting galaxies. 



3.3 Dust 



Our calculations have ignored dust. To understand the im- 
pact of dust on the Lya radiation fleld requires understand- 
ing the Lya transfer process through the ISM of galax- 
ie^j. This is because the scattering process causes Lya pho- 
tons to traverse different path s through the ISM than rest- 
fram e UV continuum photons dcharlot fc Falllll99ll : iNeufeldl 
Il99ll : [Hansen fc Ohll2006l : iLaursen et all 120091 ). As a result 



3.4 Peculiar Velocities 



dust has a different impact on Lya line and UV continuum 
photons. Observations indicate the Lya escape fraction de- 
creases with the dust content of galaxies. More specifically, 
the mean Lya equivalent width decrea ses with the observe d 
reddening of Lyman B reak Galaxies (jShaplev et al.|[200^ . 
IVerhamme et al.l (|2008l l have used their models to quantify 
the actual escape fraction of Lya photons from galaxies, 
/d.eac, as a fuuctiou of the observed dust reddening, and 



^ Several groups have included dust wh en modeling the observed 
number density of Lya emitters (e.g . |Haiman fc Spaanj [l999l : 
iMao et al.ll2007l : iKobavashi et al.ll201(]l : IPaval et al.ll2010r i. When 

comparing such models with actual data, a degeneracy exists be- 
tween the opacity of the ionized IGM and dust to Lya radiation 
(i.e. predicted number densities are not sensitive to what the ex- 
act source of opacity is). As argued in this section, understanding 
how outflows shape the intrinsic Lya emission line that emerges 
from galaxies provide such models with important new constraints 
on the IGM-and therefore-dust opacity. 



10 galaxy by 


Pello et al.l ll2004l)-which was 


Bremer et al. 


20041: Weatherlev et al.l |200 



ICen et al.l 



neutral IGM, if the Lya emitting region in the galaxy was 
receding relative to the surrounding absorbing IGM with 
V ^35 km s~^. In this study, the galaxy ionized an H II re- 
gion with a radius of ~ 0.28 pMpc. While peculiar velocities 
of this magnitude do not play an important role in reducing 
the Gunn-Peterson damping wing optical depth (for which 
one needs Av J>500 km s^^, Eq[2|, they can play a role in 
reducing the opacity of the ionized 'local' IGM (see § 13. ip . 
just like outfiows. However, both mechanisms are expected 
to leave different signatures on the observed Lya line shape: 
scattering off HI outflows is expected to result in asym- 
metric emission lines (see Fig[2|, while peculiar velocities 
ma y leave an i ntrinsically symmetric emission line symmet- 
ric (|Cen et al.|[20 05). We caution that symmetric Lya emis- 
sion lines may be observed as a result of scattering through a 
static ISM (see Fig [2]), but point out that in the latter case 
the Lya peak flux density is redshifted by ~ 500 km s~^ 
relative to the galaxy's systemic redshift. This shift is much 
larger than that expected for the model that uses peculiar 
velocities, which is ~ 35 — 70 km s~^. 
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3.5 Comments on the Model 

The model used to generate Lya line profiles in this pa- 
per assumes a single, thin, compact, sph erical shell of HI 
gas surrounding the Lyo: source (following I Ahn et al.|[20o3 : 
Verhamme et al." 2006l . l2008l : IVanzella et al.ll 



20091 ). Recently. 



Stcidcl ct al. (2010|) described an alternative simple outflow 
model that folds in constraints from observed proflles of in- 
terstellar metal absorption lines. This model can also repro- 
duce observed Lya line shapes in LBGs. In this model, the 
outflow can extend out to ~ 100 kpc, and outflow velocity 
increases with radius. It is not clear whether the models are 
incompatible: one can imagine the outflow extends over a 
wide range of velocities and spatial scales, but that the out- 
flow's opacity to Lya photons is dominated by gas closer to 
the galaxy that is restricted to a narrower range of veloci- 
ties. As a result, it is possible to model the impact of the 
outflow on the Lya radiation held as a single shell that has 
a low velocity dispersion. This is clearly an issue that needs 
further investigation. However, a more detailed investigation 
is beyond the scope of the present paper, especially because 
our main results are not likely sensitive to the precise model 
that one uses to describe the impact of the outflow on the 
Lya radiation fleld. 

In detail, the directly transmitted fraction of Lya pho- 
tons through the IGM at the highest redshifts (the 'flrst' 
galaxies a.t z > 10) probably depends quantitatively on the 
precise outflow model that one uses, since different mod- 
els may imply different predicted redshift evolution. For ex- 
ample, in the models of I Verhamme et al.l (2006, 200^ ) and 
IVanzella et al.l l|2009l ) the emission in the red wing of the 
line (at Av ^500 km s~^) comes from Lya photons that 
scatter back and forth repea tedly between the e xpanding 
HI outflow. In the model of ISteidel et al.l (|2010l ') this fre- 
quency 'diffusion' does not occur, and the emission in the 
red wing of the line reflects the maximum velocity in the 
outflow. In this latter model, the flux in the red wing of 
the line is suppressed at z > 10 if the outflow velocities 
are much lower at high redshift (see the end of § 12. 3[) . On 
the other hand, it is difficult to completely suppress fre- 
quency diffusion. Significant frequency diffusion occurs even 
in clumpy outfiows w ith a covering factor of only fcov ~ 0.5 
(jHansen fc OhI l2006l . their Fig 19), which is less than the 
outflow covering factor inferred bv ISteidel et al ] (|2010l ) over 
a signiflcant range of radii. We therefore expect that photons 
likely emerge at Aw ^500 km s~^ from the 'flrst' galaxies in 
both models, but that the fraction of the flux that emerges 
at these frequencies-and therefore /trans-depends on the de- 
tailed properties of the HI in outflows in galaxies at these 
redshifts. This underlines one of our conclusions presented 
in § |4] (also see § 13. 1|) . that it is crucial to understand the 
gas kinematics in the ISM (and/or circumgalactic medium) 
of high-redshift galaxies in order to infer the properties of 
the IGM from the observed Lya flux. 

Our main results regarding galaxies during the later 
stages of-and after-reionization {z <.8) require less (or no) 
extrapolation to higher redshifts, and practically do not de- 
pend at all on the precise underlying model of the Lya line 
shape. Instead these results rely only on the fact that the 
observed Lya line shape can extend to freq uencies that lie 
well b eyond Av ~ 500 km s~^ out to z = 5.5 jVanzella et al.l 
I2OO9I ). 



It has been argued that radiation pressure on dust 
grains-which in turn are coupled to the interstellar gas- pro- 
vide an important source of pressure in the ISM of galax- 
ies, and that the observed outflows in galaxies are driven 
predominantly by radiation pressure (see e.g. Murray et al. 
2005). Given that the highest redshift galaxies likely con- 
tained little dust (see 13. 3[) . radiation pressure may be less 
important in these galaxies, and outflows could be weakened 
signiflcantly. Interestingly, especially in these flrst galaxies, 
a signiflcant fraction (~ 20%) of the total bolometric lumi- 
nosity of the galaxy is Lya line radiation. If this radiation 
is 'trapped' by a large column of H I gas, then the radiation 
pressure exerted by this Lya radi ation itself becomes im - 
portant in driving the H I gas out (|Diikstra fc Loebll2008l '). 
thereby enhancing the detectability of the Lya emission. 



4 CONCLUSIONS 

The next generation of telescopes aim to directly observe the 
flrst generation of galaxies that initiated the reionization of 
our Universe. The Lya emission line is robustly predicted 
to be the most prominent intrinsic spectral feature of these 
galaxies. In this paper we investigated the prospects for de- 
tecting this Lya emission, taking account of radiative trans- 
fer effects that are likely to occur in the interstellar medium 
(ISM) of these galaxies. 

Observed interstellar metal absorption lines (Si II, O I, 
C II, Fe II and Al II) in Lyman Break Galaxies (LBGs) 
are typically strongly redshifted relative to the galaxies' sys- 
temic velocity (with a median off-set of ~ 160 km s~^), while 
the Lya emission line is strongly redshifted (with a median 
velocity offset of ~ 450 km s~^ iSteidel et al.ii201Q . also see 
Shapley et al. 2003). Thi s suggest that large scale outflows 
are u biquitous in LBGs jShaplev et al.l l2003l : ISteidel et al.l 
I2OIOI ). Furthermore, scattering of Lya photons by H I in out- 
flows has successfully explained the observe d Lya line shapes 
in Lya emitting galaxies at z = 3 — 6 (e.g. IVerhamme et all 
l2006l . I2OO8I : IVanzefla et al]|2009l ). Scattering off outflows of 
interstellar H I gas can shift Lya photons to the red side 
of the line before it reaches the IGM (Fig [3]). At these fre- 
quencies the Gunn-Peterson optical depth may be reduced 
to order unity as a result. 

In this paper we investigated the detectability of Lya 
radiation under the assumption that the outflows observed 
at low redshift also occur in the highest redshift galaxies. We 
found that outflows may cause as much as ^5% of the emit- 
ted Lya radiation to be transmitted directly to the observer, 
even through a fully neutral IGM (ij |2.3[) . Since the intrinsic 
(restframe) equivalent width of the Lya line can be as high 
as EWint — 1500 Afor the flrst generation of galaxies, the ob- 
served EW= /transEWint = 75(/tran. /0.05) (EWint / 1500 A) 

A. For compariso n, only 4% of the z = 3 LBG population 
have larger EWs (|Shaplev et all 120031 '). We showed that for 
/trans J>3% it may be easier to detect Lya line emission with 
NIRSPEC on the James Webb Space Telescope (JWST), 
than continuum radiation with JWST's NIRCAM in the 
same integration time. We also note that the next genera- 
tion of ground-based 30-m telescopes with diffraction limited 
AO are expected to be more (less) sensitive than JWST at 
A ^1.3/im for R 3> 100 (R <^W0) spectroscopy (see Moun- 
tain et al. 2009, their Fig 2). That is, ground based high- 
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Figure 6. Schematic explanation for wliy outflows promote tlie detectability of Lya emission from galaxies surrounded by signiflcant 
amounts of neutral intergalactic gas: In the top panel a galaxy is surrounded by a large bubble of ionized gas, which in turn is surrounded 
by neutral intergalactic gas. Lyo emission from this galaxy redshifts away from resonance as it propagates freely through the H II bubble 
(as indicated by the line color). Once the Lya photons reach the neutral IGM, they have redshifted far from resonance where the Gunn- 
Peterson optical depth is reduced tremendously (see Eq[2]l. Because of the reduced CP optical depth, some fraction of the emitted Lya 
is transmitted to the observer without scattering in the IGM. In this drawing, the thickness of the line represents the specific intensity 
of the Lya radiation field. The bottom panel shows that outflows surrounding star forming regions (represented by the expanding ring. 
The far side is receding form the observer and has a larger redshift, which is represented by the color) can Doppler boost Lya photons 
to frequencies redward of the Lya resonance. In this scenario, a fraction of Lyo can propagate directly to the observer without the H II 
bubble. 



resolution spectroscopic searches for high redshift galaxies 
can detect fainter galaxies at a fixed integration time when 
/trans J>5%. Such scarches can therefore be competitive with 
searches that employ the drop-out technique. Irrespective of 
the survey strategy that is used to search for the highest 
redshift galaxies, the prospect that Lya can provide galax- 
ies with spectroscopic redshifts is promising, and important 
as no robust predictions exists for the detectability of other 
emission lines (fj 13. 2p . 

This paper has focused on the first generation of galax- 
ies that were surrounded by a neutral IGM, but our work 
also applies more broadly. For example we argued in § 13.11 
that H I outflows promote the detectability of the Lya emis- 
sion line during later stages of reionization when much of the 
absorption is resonant absorption in a highly ionized HII re- 
gion. As a result outflows can reduce the minimum H II 
bubble size that is required to render LAEs 'visible'. This is 
illustrated schematically in Figure [6l Similarly, H I outflows 
also promote the detectability of the Lya emission line after 
reionization has been completecjf] I3.1|) . 

° Recently, IZheng et al. I ll2010d) showed that scattering in the 



In summary, radiative transfer effects in the ISM of high 
redshift galaxies have been shown to broaden, and-in the 
case of outfiows-redshift the emergent Lya fiux to a level 
that allows 5% or more of the photons to escape absorption 
by the neutral IGM. Goupled with the large intrinsic Lya 
line EW of the first generation of galaxies, we have shown 
that searches for galaxies in their redshifted Lya emission 
line can be competitive with the drop-out technique out to 
the highest redshifts that can be probed observationally in 
the JWST era. 
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'local' IGM immediately surrounding Lya sources can introduce 
a unique anis otropy in the t wo-poi nt correlation function of LAEs 
at z = 5.7. IZheng et all l l2010hl) argued that this scattering- 
induced signature is reduced when the intrinsic (i.e. prior to scat- 
tering) Lya line width is enhanced. Winds are therefore also ex- 
pected to reduce this clustering signature. That is, the clustering 
of LAEs post-reionization can provide constraints on the impor- 
tance of winds in shaping the Lya line shape. 
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Figure Al. Left panel: same as Figure [S] Tiie red histogram siiows tiie spectrum of thie 'point source' (unscattered radiation) as 
extracted from the Monte-Carlo code. Right panel: the solid black line shows the surface brightness profile of the standard Loeb-Rybicki 
halo (scattered radiation, as in Fig[TJ, while the red dashed line shows the surface brightness profile of the scattered radiation for the 
outflow model. In the outflow model, the surface brightess profile of the Lya halo is boosted significantly at 9 <^5 arcsec. However, this 
boost goes away when a small (radius=50 pkpc) is present around the galaxy (as shown by the blue dotted line). This figure shows that 
while outflows boost the dctectability of the point source, they do not boost the detectability of the Lyo LR-halos. 

APPENDIX A: LOEB-RYBICKI HALOS WITH OUTFLOWS 

In the main body of this paper we focused on Lya radiation that was transmitted directly to the observer through a fully 
neutral IGM. This radiation is confined to an angular region that is set by the physical scale of the outflow (throughout the 
paper we referred to this as the 'point source'), and is much more easily detectable than the scattered Lya radiation that is 
in a diffuse halo. For completeness we show how the surface brightness profile of the scattered radiation is modified because 
of the outflow in the galaxy. 

In the left panel of Figure lAll the black solid line shows the spectrum of unscattered radiation for the model with 
(A'hi, Vsh) = (10'^° cm~^,200 km s^^) (as in Fig|3l The grey solid line is again the intrinsic spectrum that emerges from the 
galaxy). The red histogram shows the spectrum of the point source as extracted from the Monte-Carlo code. The agreement 
between the Monte-Carlo and analytic solution (given by Eq|4]) are again good. 

In the right panel the black solid line shows the surface brightness profile of the standard Loeb-Rybicki halo (as in Fig[TJ, 
while the red dashed line shows the surface brightness profile of the scattered radiation for the outflow model. Clearly, the 
surface brightess profile is boosted signiflcantly at <.5 arcsec. Radiation that escaped from the galaxy far in the red wing 
of the line (at Av ^500 km s^^) is most likely to scatter in close proximity to the galaxy, where it is closest to resonance. 
After this scattering event, these photons again have a non-negligible probability of propagating directly to the observer. This 
boosts the surface brightness in the inner region of the Lya halo. We point out that the maximum boost of a factor of ~ 10 
is reached at 9 <.l arcsec. In reality we expect gas to be ionized at such short distances from the galaxy. Indeed, if we insert 
a small (radius = 50 pkpc) H II bubble around the galaxy, the we find that the boost disappears, and we almost recover 
the original surface brightness profile (as indicated by the blue dotted line). Larger H II regions only suppress the surface 
brightness proflle of the central core. The main point of this calculation is to demonstrate that the radiative transfer processes 
that we invoked to boost the detectability of the 'point source' do not also boost the detectability of the Lya halos. 
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